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Abstract 

Canonical ensemble molecular dynamics simulations of liquid methanol, modeled using a rigid- 
body, pair-additive potential, are used to compute static distributions and temporal correlations of 
tagged molecule potential energies as a means of characterising the liquid state dynamics. The static 
distribution of tagged molecule potential energies shows a clear multimodal structure with three distinct 
peaks, similar to those observed previously in water and liquid silica. The multimodality is shown to 
originate from electrostatic effects, but not from local, hydrogen-bond interactions. An interesting 
outcome of this study is the remarkable similarity in the tagged potential energy power spectra of 
methanol, water and silica, despite the differences in the underlying interactions and the dimensionality 
of the network. All three liquids show a distinct multiple time scale (MTS) regime with a 1//" 
dependence with a clear positive correlation between the scaling exponent a and the diffusivity. The 
low-frequency limit of the MTS regime is determined by the frequency of crossover to white noise 
behaviour which occurs at approximately 0.1 cm~^ in the case of methanol under standard temperature 
and pressure conditions. The power spectral regime above 200 cm~^ in all three systems is dominated 
by resonances due to localised vibrations, such as librations. The correlation between a and the 
diffusivity in all three liquids appears to be related to the strength of the coupling between the localised 
motions and the larger length/time-scale network reorganizations. Thus the time scales associated with 
network reorganization dynamics appear to be qualitatively similar in these systems, despite the fact 
that water and silica both display diffusional anomalies but methanol does not. 
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charusQchemistry . iitd . ernet . in) 
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1 Introduction 



Considerable experimental as well as theoretical attention has been devoted to characterizing the 
structure and dynamics of hydrogen-bonded liquids [ll[2l[3llll[5l[6l[7l[8l[9l [10] . The motivation for 
these studies is primarily due to the importance of water as a solvent for biological and chemical 
processes. A comparison of the behaviour of water with other second-row hydrides, such as HF, NH3 
and CH3OH, is essential for understanding the nature of hydrogen-bonded dynamics and the extent to 
which the behaviour of water is unique. Methanol, in particular, is of interest because it is the simplest 
molecule which can exhibit both hydrogen bonding and nonpolar interactions. 

While interactions in simple liquids are dominated by steep, short-range repulsions and long-range, 
isotropic attractions, hydrogen-bonded liquids have strong, local anisotropic interactions. In the case 
of water, each molecule can form at most four hydrogen bonds, leading to a three-dimensional, open, 
locally tetrahedral network structure. The strength of hydrogen bonds is estimated to lie between 5 
and lOksT at melting [llj which is strong enough that a substantial fraction of hydrogen bonds will 
be intact at room temperature. Thermal fluctuations will, however, be large enough, in comparison 
to the bond strength, to ensure that such bonds will have a finite lifetime of the order of picoseconds. 
As a result, the dynamics of the liquid will be dominated by the behaviour of the three-dimensional, 
hydrogen-bonded network, parts of which are constantly broken and reformed. In the case of methanol, 
hydrogen bond strengths are similar to those in water [12j, but each methanol molecule can form at 
most three hydrogen bonds, of which only one can be as a proton donor. In liquid methanol, however, 
simulations as well as neutron scattering experiments show that the typical number of hydrogen bonds 
per molecule is two \13\ [T^ 115 ]. Consequently, linear chains with very rare branch points are seen, 
rather than the three-dimensional network characteristic of water. 

Hydrogen-bonded liquids can be thought of as a subset of network-forming liquids, other examples 
being ionic melts such as ZnCl2, Si02 and BeF2 [16j. The strong local coupling of individual atomic or 
molecular units leads to the existence of multiple length and time scales corresponding to cooperative 
rearrangements of the network involving different numbers of molecules. To characterise the dynamical 
behaviour of network-forming liquids, it is convenient to use the power spectral density of a mechanical 
quantity that is sensitive to motion on a number of different length scales. The power spectral density 
of an observable A(t) as a function of time t over a time interval T, is defined as|17] 



T^oo T J-T/2 

where S{f) is the spectral power associated with frequency /. In a recent set of studies, we have 
shown that the power spectrum associated with the tagged particle potential energy can be used 
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to characterize the dynamics of tetrahedral, network-forming hquids, such as water [T8| \T9\ [20l [2H 
[22] and Si02 [23J. Simple hquids show a white noise power spectrum where S(f) is independent of 
frequency. Network-forming liquids, however, have a characteristic multiple time-scale regime with 
1//" dependence on frequency / with a lying between 1 and 1.8 [2jj. The three liquids that we have 
studied all show a region of anomalous diffusivity where mobility increases on isothermal compression; 
this density-dependent variation in mobility is strongly correlated with the scaling exponent a. Our 
results therefore suggest that monitoring the power spectrum of tagged particle potential fluctuations 
to be a simple and direct method for linking phenomena on three distinctive length and time scales: 
the local molecular environment, hydrogen bond network reorganisations and the diffusivity. 

In this work, we characterise the multiple time-scale behaviour of liquid methanol, focusing primar- 
ily on the power spectrum of the tagged particle potential energy. Methanol presents an interesting 
contrast to the systems we have studied earlier since hydrogen bonding results in formation of one- 
dimensional chains, rather than a three-dimensional network. Unlike water, methanol does not show 
any anomalous liquid state properties, such as the diffusional anomaly mentioned above, or ther- 
modynamic anomalies related to the existence of a temperature of maximum density under isobaric 
conditions. It is therefore interesting to compare the manner in which the multiple time-scale be- 
haviour of methanol and water differ. We perform molecular dynamics simulations of methanol using 
Haughney's HI potential |25j which predicts thermodynamic, structural and transport properties in 
fair agreement with a wide variety of experimental data |131 [T^ [26l [271 [28] as well as recent ab initio 
molecular dynamics simulations |29[ [30j . Section 2 contains details of the molecular dynamics simu- 
lations as well as of the power spectral analysis. Section 3 contains a discussion of our results, and 
Section 4 presents the conclusions. 

2 Computational Details 

2.1 Potential energy surface 

The intermolecular interactions in methanol are modeled using Haughney's HI pair-potential which 
treats methanol as a rigid three-site entity where the centres of interaction are positioned on the carbon 
(C), oxygen (O) and the hydroxyl proton Ho[25j. The interaction between any two methanol molecules 
a and /3 has the following functional form: 



where the indices i and j label the atomic sites on the two molecules and corresponds to the site-site 
separation. The atomic sites are labelled by the partial charge, and the Lennard-Jones parameters an 
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and en, which are summarised in Table I. The Lennard- Jones parameters for interactions between unhke 
atoms are obtained using the Lorentz-Berthelot mixing rules: aij = {ai + cTj)/2 and eij = ^e-i^y The 
molecular geometry is specified using the following bond lengths and bond angles: r(C-0)=1.4246^, 
r(0-Ho)= 0.945li, r(C-H)=1.0936i, ZCOHo= 108.53° and ZHCH =108.63°. The methyl protons 
are not considered as interaction sites but they are included when specifying the molecular geometry, 
so that the molecules have correct masses and moments of inertia. 

2.2 Molecular Dynamics 

Molecular Dynamics (MD) simulations of a system of 256 methanol molecules were carried out in canon- 
ical (NVT) ensemble, using the DL_POLY software package [3T1[32], under cubic periodic boundary 
conditions. The effects of electrostatic, long-range interactions were accounted for by the Ewald sum- 
mation method. The non-coulombic part was truncated at half the box length. Berendsen thermostat, 
with a time constant rB=200ps, was used to maintain the desired temperature for the production run. 
The Verlet algorithm with a time step of 2 fs was used to integrate the equations of motion. The 
quaternion algorithm was used to maintain the rigid body constraints on the methanol molecules. 

The system was studied along two isochores, 0.768 g cm~^ and 0.878 g cm~^. The experimental 
density at 300 K and latm is 0.768 g cm~^ [25] while 0.878 g cm^'^ is the density around which the 
glass transition was reported to occur at around 155 K using the HI potential [33]. The temperatures 
studied were 300, 250, 210, 180, 170 and 155 K for the 0.768 g cm^^ isochore and 300, 250, 210, 
180 and 170 K for the 0.878 g cm^^ isochore. The initial configuration for the molecular dynamics 
simulations at 300 K was generated using the PACKMOL program [34J. The system was then slowly 
cooled to obtain lower temperatures. The system was equilibrated for 1-4 ns for the temperature range 
of 300-155K, and the production run lengths were kept between 2-8 ns. The details of run lengths 
are summarised in Table II. Self-diffusivities were computed using the Einstein relation [35j. Long 
run lengths were used, particularly at low temperatures, to ensure adequate sampling of configuration 
space, as measured by the mean squared displacement (MSD) of oxygen atoms, as shown in Figure 1. 
The time of onset of the diffusional regime, id*//) was identified from the MSD versus time plots as the 
regime in which the Einstein relation is obeyed. Along an isochore, as the temperature is decreased we 
observe that the onset of diffusional regime shifts to higher and higher times (0.5 to 2.5 ns). At 155 K, 
for the 0.768 g cm^'^ isochore, no linear dependence of the MSD on time is observed even upto 8 ns run 
length, suggesting that the system is out of equilibrium and we are dealing with an arrested structure. 
Table II shows the time of onset of diffusional regime, ^di//, as determined from the applicability of 
the Einstein relation to describe the time dependence of the MSDs. Our simulation results are in 
agreement with previous results [25^ [33] within the statistical error bars. 
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2.3 Generating Power Spectra 

Based on our previous work on water where we considered various tagged particle observables, such 
as local orientational order metrics and tagged particle potential and kinetic energies, we concluded 
that the tagged particle potential energy is sensitive to a wide range of time scales associated with 
the hydrogen-bonded network [19j. Therefore, in this study, we focus on the static distributions 
and power spectra associated with tagged molecule potential energies, Su{f)- The tagged molecule 
potential energy, u{t), corresponds to the interaction energy of an individual molecule with all the 
other molecules in the system. Since the configurational potential energy function is pair-additive, the 
total potential energy, U{t) = 0.5 J2iUi{t) where the sum extends over all molecules. 

The tagged particle potential energies were stored at intervals of 10 fs i.e. every 10 MD steps, 
which corresponds to a Nyquist frequency of 1666 cm~^. The value of the Berendsen thermostat time 
constant (r^) provides the lower limit on the frequency range over which we can obtain reliable power 
spectra; thus, tb=200 ps corresponds to a lower frequency limit of 0.165 cm^^. Standard fast Fourier 
transform routines were use with a square sampling window. The normalization convention was chosen 
such that the integrated area under the S{f) curve equalled the mean square amplitude of the time 
signal. Windows containing 2^^ data points were used for Fourier transformation [36j. Statistical noise 
in the power spectra was reduced by averaging overlapping time signal windows as well as individual 
tagged molecule spectra [T5] . 

Our analysis of the power spectra associated with tagged molecule potential energy focuses on 
the three distinct regions of the spectrum: (i) identification of resonances or power spectral peaks, 
typically corresponding to high-frequency local vibrational modes which are weakly coupled to network 
reorganizations; (ii) identifying the frequency range and the exponent associated with the l/f°' multiple 
time scale regime, and (iii) the frequency of crossover to white noise behaviour. A linear least-square 
fitting of ln5(/) was used to obtain the scaling exponent a associated with the multiple time-scale 
region. 

2.4 Model-based fitting procedure for power-law-like spectra 

A procedure for quantitatively fitting the entire power spectrum, and not just the multiple time scale 
regime, to a physically motivated parametric model for the time-scale distribution associated with 
molecular motions in water has been developed previously [22]. In this paper, we apply a similar 
procedure to quantify the temperature-dependent changes that take place along the 0.768 and 0.878 
g cm~'^ isochore of HI model of liquid methanol. The analytical model is a sum of six functions that 
correspond to well-defined model behaviors: 
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• a Sp{f]a,j3) = a? / f"^ power-law noise term containing two parameters, corresponding to the 
amplitude and spectral index; 

• five resonance distributions. Here we have taken a fiat distribution of resonance frequencies lying 
between fmin and fmax with fixed relaxation rate A. The resulting spectral densities are 

/fmax r „ „ „1 — 1 

A -|-47r (/ - /o) J dfo. (3) 
J min 

For each resonance, four parameters, corresponding to the amplitude, minimum and maximum 
resonance frequencies and relaxation rate, must be introduced. 

In addition the model function also contains a white noise term and a filter function 

^^^'''"^= ((2.//aV + 1) 
with n 5. The filter function takes into account the steep fall off in the spectral density at high 
frequency, presumably due to the upper limit imposed by the time step of the molecular dynamics 

algorithm. The filter function is approximated rather roughly by the function F(f;X,5), the numerical 
tests show that that the fit is rather insensitive to its exact shape. The overall 24-parameter model is 
given by: 



Smodelif) = F{f;Xo,5) 
where C is an additional white noise term. 



fc=6 

Spif'i Ol, + Sresif'i O-k-i fmin,k-i fmax,ki Afc) 
k=2 



+ C (5) 



3 Results and Discussion 

3.1 Static distribution of tagged molecule potential energy 

Figure 2 shows the static distribution, P('u), of the tagged molecule potential energies, u{t), for 
methanol at several temperatures along the 0.768 and 0.878 g cm~^ isochores. Unlike in the case 
of monoatomic van der Waals liquids, methanol resembles water and silica in that, it shows multi- 
mododal P{u) distributions indicating that there are energetically distinct local environments in the 
liquid. The presence of distinct peaks suggests that there is a small energetic barrier to transit from one 
energetic environment to another. This heterogeneity in local environments seems to be characteristic 
of network-forming liquids, regardless of the dimensionality of the networks. It should be noted that 
the number of peaks in the P{u) distributions of water and silica are larger but the individual peaks 
are less distinct. 

The three peaks in methanol are centred at approximately at ~ -130, -85 and -30 kJ mol~^. 
The distinct demarcation between the three peaks as well as the relative probabilities change with 
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temperature, with the lowest energy peak becoming less prominent. At low temperatures, just above 
the glass transition, the relative intensities of the peaks are very sensitive to the quality of equilibriation. 
The mean value of the tagged potential energy, (u), varies between -82.4 and -69.2 kJ mol~^ between 
155 and 300 K along the 0.768 g cm^^ isochore. A similar behaviour is observed in case of 0.878 g 
cm~^ isochore. The mean value of the tagged potential energy, (u), ranges from -84.17 to -72.3 kJ 
mol~^ between 170 and 300 K along this isochore. 

Figure 3 shows the probability distributions corresponding to the tagged molecule potential energy 
contributions from the van der Waals, screened Coulombic and reciprocal space contributions to the 
total potential energy for a single state point (300K, 0.768 g cm^^). It is immediately obvious that the 
multimodal character is due to the electrostatic contribution from the screened Coulombic interactions. 
The P{urec) and P[uyd\^) distributions are unimodal and very narrow in comparison with -^('U^creen)* 

The dependence of the multimodality of the static tagged particle potential energy distribution 
on the electrostatic contribution was also noted in our previous work on silica [23J. In the Haughney 
model, the Coulombic contribution is parameterised to take into account the net effect of both hydrogen 
bonding and multipole moment interactions. To understand the extent to which the peak structure 
of the P{u) distributions is correlated with the hydrogen bonding, we have classified the methanol 
molecules based on their tagged particle potential energy into three categories, as defined in Table 
III, and computed the fraction /„ of n-bonded hydrogen molecules in each category. The geometric 
criteria for identifying the hydrogen bonds were taken from ref.[25]. Two molecules are regarded to be 
hydrogen bonded if roo < 3.5^, tohq ^ 2.6^ and the angle IHqOO < 30°. The simulation averages 
of /o, fi, f2 and /„ at 300K and 0.768 g, shown in Table III, reproduce those of Haughney et. al. 
within statistical error. Activation energes of hydrogen bond lifetimes was found to be 15.6 kJ mol~^, 
which is close to the previously reported values of 16 kJ mol~^ \25\ I33j . 

The results in Table III show that non-bonded molecules are essentially absent in the 170K to 300K 
regime. The peak centered at -30 kJ mol~^ has a somewhat higher value of /i and the peak centered 
at -130 kJ mol^^ has a somewhat higher value of /s than the overall simulation averaged values of /i 
and /a- The CH3OII molecules which belong to the central peak of the P{u) distribution have an /„ 
distribution which is essentially identical to the overall ensemble averaged distribution. Therefore the 
peak structure of the P{u) distributions cannot be attributed to differing degrees of hydrogen bonding, 
though there is a weak correlation between Uhi, and u. Our previous work on water, as well as the 
results in Section 3.4, show that as the electrostatic contribution is turned off and the Lennard-Jones 
component relatively enhanced, the degree of multimodality diminishes. The multimodality must 
therefore originate from longer range electrostatic interactions, stretching over at least next nearest 
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neighbours, in addition to nearest neigbours. It would be interesting, in this context, to study fluids 
with dipolar, and perhaps higher order multipolar, interactions [35j . 

3.2 Power spectra of tagged particle potential energy 
3.2.1 Comparison of methanol and water 

Figure 4 compares the power spectra generated by the fluctuations in the tagged potential energy of 
methanol, £'„(/), at 300 K and 0.768 g cm^^ ^i^^ ^^^^^^ SPC/E water at 1.0 g cm'^ at 300 K. The 
densities correspond to those observed under ambient temperature and pressure conditions. Despite 
the difference in the dimensionality of the hydrogen bonded network between methanol and water, 
the key power spectral features show strikingly similar qualitative features with: (i) high frequency 
localised modes (ii) a 1//° multiple time scale regime with exponent a at intermediate frequencies and 
(iii) a crossover to white noise at low frequencies. The high-frequency cutoff of the power spectrum is 
determined by the time step used in the molecular dynamics simulations. 

Given the rigid molecule approximation used in the present simulations, no intramolecular vi- 
brational modes can be observed in our studies. We briefly summarise the information on speciflc 
intermolecular vibrational modes available from previous experimental and simulation studies. In the 
case of water, there is a broad librational band between 500 and 800 cm~^, as well as two- molecule 
0-0 stretching and three-molecule 0-0-0 bends at approximately 200 cm^^ and 60 cm~^ [37 1 1381 [39] , 
In the case of methanol, experimental IR spectra locate peaks at approximately 680, 270, 135 and 65 
cm~^ \28 \ 1401 HT] . The hindered rotations of the methanol molecule, with the methyl group treated as 
a united atom in the OPLS model [12], are associated with frequencies of approximately 580, 210 and 
50 cm~^ f32]. Explicit treatment of the geometry of the methyl group in the HI potential shows that 
the 270 cm~^ peak is associated with the velocity autocorrelation function of the hydroxyl hydrogen 
atom [25]. More recent studies of the optical spectra of methanol identify the COH pure librational 
peak as lying between 200 and 300 cm^^ and the peak at approximately 600 cm~^ as the H-bond 
constrained libration around the CO bond [431l44j. The low frequency librational bands at 65 and 135 
cm~^ are expected to be strongly coupled with intermolecular vibrational modes. 

The Su{f) spectrum of water shows a broad librational band between 500 and 800 cm^^ but 
the intermolecular stretches are not observable at this density and temperature since they lie in the 
multiple time-scale (MTS) region. In the case of methanol, the two high frequency librational bands 
are observed as distinct peaks at 270 and 600 cm~^. Lower frequency vibrational bands at 65 and 135 
cm~^ are not visible as distinct peaks in the Su{f) spectrum of CII3OH at this state point because 
they lie in the MTS region. 
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The multiple-time scale (MTS) regime in methanol at 300 K and 0.768 g cm ^ extends from 0.2 to 
200 cm~^. There is, however, a small but distinct change in the exponent a at 40 cm~^; the a value 
for 0.2-40 cm-i range is 1.36 while for 40-200 cm'^ it is 1.22. For SPC/E water at 300 K and 1.0 g 
cm~^, the 1//" region extends from 1-200 cm^^ with a = 1.4. The presence of a larger number of 
high-frequency resonances in methanol, compared to water, maybe responsible for this change in the 
exponent of the MTS regime at approximately 40 cm^^. In the remainder of this paper, we refer as 
the MTS regime to the frequency interval lying below 40 cm^ and above the frequency of crossover to 
white noise. 

At the equilibrium densities associated with 1 atm pressure and 300K, the white noise regime in 
water occurs below 1 cm~^ whereas for methanol, it begins at a lower frequency of 0.2 cm~^. This 
suggests that the time taken for the potential energy fluctuations to decorrelate is smaller for water 
than methanol and therefore the overall relaxation dynamics over length scales of a few molecules will 
be faster in water than in methanol. This is consistent with the viewpoint expressed by Ladanyi and 
Skaf [1], based on orientational correlation functions, angular VACFs and NMR relaxation rates, that 
water and methanol relax at similar scales at short times but the relaxation of CH3OH becomes much 
slower for times greater than 0.6ps (55.5 cm~^). The overall faster dynamics of water, despite the 
existence of a stronger hydrogen bond and an extended three-dimensional network is attributed to the 
higher translational diffusity due to lighter mass, more rapid librational modes and the coincidence of 
the direction of the molecular dipole moment with one of the principal axes of rotation [1] . 

The power spectra associated with the fluctuations in the total tagged potential energy and the 
contributions from reciprocal and screened energies arising out of Ewald summation and the van der 
Waals interactions were computed at 300K and 0.768 g cm~'^. Though the figure is not shown, the 
behaviour of is almost exactly reproduced over much of the frequency range by the screened 

charge contribution. The spectral power associated with the reciprocal space and van der Waals 
contributions is lower and, therefore has a minor effect on the overall spectrum except in the high- 
frequency region. The u^dw contribution produces a relatively unstructured power spectrum, more like 
Argon, while the Coulombic interactions are responsible for much of the high-frequency, short-length 
scale structure. This is consistent with our earlier observations on silica |23j . 

3.2.2 Temperature and Density Dependence 

The temperature dependence of the power spectra associated with tagged molecule potential energies 
is considered in Figure 5. Figures 5(a) and 5(b) show 5„(/) curves for different temperatures along the 
0.768 and 0.878 g cm~'^ isochores respectively. The glass transition, as reflected in the time-dependence 
of the mean-square displacement plots, occurs at a lower temperature for the higher density isochore. 
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The Su{f) plots for both densities at 300 K show well-defined librational peaks. As discussed in 
section 3.2.1, the 1//" MTS region can be partitioned into a low-frequency MTS regime extending 
from 0.2-40 cm~^ with a„ = 1.36 and a higher frequency MTS regime from 40-200 cm~^ with exponent 
1.22. 

At this temperature, there are no peaks lying between 40 and 200 cm~^ though molecular centre- 
of-mass VACFs, far-IR and dielectric loss spectra suggest that vibrational modes at 65 and 135 cm^^ 
corresponding to intermolecular vibrational modes involving two or three molecules should be present. 
The change in the value of a at approximately 40 cm~^ suggests that these two- or three- molecule 
network modes are partially decoupled from the longer length scale, lower frequency network reaar- 
rangements. This discontinuity becomes more pronounced with the lowering of temperature and dis- 
tinct peaks at approximately 65 and 140 cm^^ can be seen in Su{f) plots at 170 and 155K. These 
peaks must correspond to caging vibrations of methanol molecules as the system dynamics slows down. 
For these low temperature isotherms, the onset of difFusional dynamics in the MSD plots is very slow 
indicating that the system is very close to glass transition. The crossover to white noise which is visible 
below 0.2 cm~^ for the 300K isotherm moves to frequencies below 0.1 cm^^ for temperatures of 250K 
or less. In this study, we did not attempt to generate power spectra at lower frequency since this would 
have implied using very large values for the time constant of the Berendsen thermostat and therefore 
would have compromised the temperature control. 

3.3 Parametric fitting of the Power Spectra 

The results of the fitting of the power spectra along the 0.768 g cm^^ isochore is shown in Figure 6. 
The quality of the fits, given the complex shape of the power spectrum, is good. It is interesting to 
compare the parametric distribution functions used for water |22j and methanol. Fitting the water 
power spectra requires only two resonances [22] while fitting the methanol power spectrum requires five 
resonance terms. In the case of both water and methanol, power spectra were computed over the 0.1 
to 1666 cm~^ regime. The low frequency white noise is more pronounced for water in this frequency 
regime than for methanol for the temperatures and densities studied here. As a consequence, a Debye 
term corresponding to fixed frequency relaxation processes is required for the water spectral fits. In 
the case of methanol, an additional white noise term and a filter function is used in the fits but no 
Debye relation term was required. 

3.4 Modifying the Strength of Hydrogen Bonding 

The Hl-model of methanol effectively condenses the information about the short-range repulsion and 
dispersion using the Lennard-Jones terms. The long- range electrostatic interactions and anisotropic 
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hydrogen bonding interactions are controlled by the partial charge distribution over the atoms. To 
study the effect of varying the strength of hydrogen bonding on the power spectra and the static 
distribution of tagged particle potential energies at 300 K and 0.768 g cm~^ , we have scaled the 
Hl-model charges. Figure 7 compares the spectra of methanol with normal charges on C, Ho 

and O with the power spectra when the charges are up-scaled and down-scaled by 10%, at 300 K and 
0.768 g cm~'^. The crossover to white-noise occurs at around 0.2 cm^^ for the system with standard 
HI model charges. When the charges are increased by 10%, we see that the librational peak becomes 
more pronounced and the crossover to white noise shifts to lower frequencies (not shown in the frequecy 
range studied) owing to the fact that now the interactions are much stronger and so is the hydrogen 
bonded network resulting in a long range correlation between the network dynamics and the local 
librations in the system. An increase in the extent of multimodality in the probability distribution of 
static tagged potential energy is also observed with the increase in charges on HI model of methanol. 
On the other hand, reduction of the charges by 10%, converts a sharp librational peak occuring in 
tagged potential energy power spectrum to almost a shoulder and induces a crossover to white noise 
at frequencies below 1 cm^^. There also occurs a destructuring or loss of multimodality in P{u) 
consistent with the weakening of the hydrogen-bonded network in liquid methanol. As the charges 
are reduced and the system approaches a Lennard- Jones fluid in its behaviour, the tagged molecule 
potential energy distribution narrows significantly indicating that the heterogeneity in local molecular 
environments responsible for multiple time-scale behaviour is lost. 

3.5 Connection between multiple time scale behaviour and diffusivity 

The self-diffusivity was computed using the Einstein relation from the MSD plots for the oxygen atoms 
(see Section 2.2). Figure 8 shows the Arrhenius plot of the diffusivity. The activation energies along 
the 0.768 and 0.878 g cm^^ isochores are 10.5 and 14.5 kJ mol"^ respectively. Haughney et al report 
an activation energy of 13 kJ mol^^ for 0.768 g cm~^ isochore; the discrepancy with our values may 
be due to the fact that our run lengths are more than two orders of magnitude longer. 

In simple liquids, diffusion occurs as a combination of collisional and cage relaxation processes. An 
additional mechanism of diffusion exists in network-forming liquids due to the coupling of local high- 
frequency modes (e.g. librations in the case of H2O) with local, two-molecule network re-organizations. 
As librational modes are more strongly coupled with network-reorganizations, the exponent a rises as 
does the diffusivity, leading to a strong positive correlation between D and as discussed in our 
earlier work on water and silica [201 123] • Figure 9 demonstrates that this relationship between Ou 
and the diffusivity also holds in methanol. Clearly, the librational coupling to larger scale H-bonded 
network reorganisations is an important diffusional mechanism in CH3OH, as in water and silica. 
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despite the differences in network dimensionality. It should be noted that we have computed a„ values 
by fitting lnS'„(/) over the 1 to 40 cm~^ regime which is expected to be particularly sensitive to the 
coupling of the two- and three-molecule modes with the network regorganisations. 

In order to make a quantitative comparison of the correlation between diffusivity and the scaling 
exponent, Figure 10 plots a dimensionless diffusivity, D*, against for CH3OH, H2O and Si02 in 
the liquid phase. The dimensionless diffusivity, D* , is defined as D* = Df?-I^ j ^fk^TJm where p 
and m correspond to the number density and molecular mass of the formula units (H2O, Si02 and 
CII3OH). The correlation is very similar in all three cases, though given the differences in intermolecular 
interactions, some scatter in the correlation plot is expected. Typically, as temperatures and densities 
increase and the network structure breaks down, the a values become less sensitive to the diffusivity; 
in the limiting case of a simple liquid, a„ is zero and uncorrelated with the diffusivity. 

4 Conclusions 

In this work, we have focused on understanding the structure and dynamics of liquid methanol by 
studying the static distributions and temporal correlations associated with fluctuations in the tagged 
molecule potential energies. While the tagged molecule potential energy is not a directly accessible 
quantity, previous work on water, aqueous alkali halide solutions and silica melt have demonstrated that 
this quantity is sensitive to network reoganisations on several time scales. In the case of methanol, we 
show that analysis of this single, frequency-dependent power spectral distribution provides an overall 
view of the dynamics, including important resonances and the multiple time-scale regime associated 
with network formation. This is particularly interesting given that the resonances at 650 cm~^, 270 
cm~^, 135 cm~^ and 65""^ have otherwise been identified using very different correlation functions 
corresponding to different experimental techniques. 

The static distribution of tagged molecule potential energies shows a clear multimodal structure 
with three distinct peaks. Similar multimodal distributions were observed earlier in the case of water 
and silica melt at low temperatures and densities. For both the systems, factors which attenuate the 
network, such as temperature, density and or ionic solutes, were found to destroy the multimodal char- 
acter and lead to unimodal Gaussian-type distributions characteristic of simple liquids. In the case of 
methanol, we have further found that the correlation between number of hydrogen bonds, as identified 
by geometric criteria, and the tagged particle potential energy is very weak. The multimodality seems 
to therefore be due to electrostatic effects, but not necessarily due to local anisotropic interactions. 
The longer range of electrostatic interactions, such as dipole-dipole interactions, presumably implies 
that the organisation of second and third neighbour shells plays a crucial role in determining the tagged 
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particle potential energy distributions. In this context, it would therefore be interesting to examine 
strongly dipolar, but not necessarily hydrogen-bonded systems. 

An interesting outcome of this study is the remarkable similarity in the tagged potential energy 
power spectra of methanol, water and silica, despite the differences in the underlying interactions and 
the dimensionality of the network. All three liquids show a distinct multiple time scale (MTS) regime 
with a I//'* dependence with a clear positive correlation between the scaling exponent a and the 
diffusivity. The low-frequency limit of the MTS regime is determined by the frequency of crossover to 
white noise behaviour which occurs at approximately 0.1 cm~^ in the case of methanol under standard 
temperature and pressure conditions. The power spectral regime above 200 cm~^ in all three systems 
is dominated by resonances due to localised vibrations, such as librations. The correlation between a 
and the diffusivity in all three liquids appears to be related to the strength of the coupling between 
the localised motions and the larger length/time-scale network re-organizations. Thus the time scales 
associated with network reorganization dynamics do not appear to be qualitatively different in these 
systems, despite the fact that water and silica both display diffusional anomalies but methanol does 
not. This suggests that the equilibrium transport properties associated with the anomalies may be 
more sensitive to the way a liquids samples the overall configuration space, i.e. the excess entropy 
scaling of transport properties, than the short time dynamics. Given that time-resolved infra-red 
techniques have been used to monitor the hydrogen-bond dynamics in water, it would be interesting 
to compare the solvation dynamics in water with that in methanol and other low-molecular weight 
alcohols. 
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Table I 

Potential Parameters of Haughney's HI model 



Atom 


e(kJ/mol) 


a{A) 




C 


0.7576 


3.861 


0.297 





0.7309 


3.083 


-0.728 


Ho 






0.431 



16 



Table II 

Summary of Computational Details for MD simulations of HI methanol. The lengths of equilibra- 
tion and production runs are denoted by tgq and tprod respectively. The time of onset of the diffusional 
regime is denoted by tdiff- 



Temperature 




0.768 g cm-^ 






0.878 g cm-^ 




(K) 


teq{ns) 


tprod ('is) 


tdiff (ns) 


teq{ns) 


tprod ('is) 


tdiff (ns) 


300 


1 


2 


0.05 


1 


2 


0.05 


250 


1 


2 


0.08 


1 


2 


0.10 


210 


2 


4 


0.20 


2 


4 


0.40 


180 


3 


6 


2.00 


3 


6 


3.00 


170 


3 


6 


2.50 


3 


6 


3.00 
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Table III 

Hydrogen bond statistics of methanol molecules at two different state points. The fraction of 
molecules with n hydrogen-bonded neighbours is denoted by /„, The ensemble averaged values of /„ 
are coomputed, in addition to the averages over molecules belong to restricted ranges in the tagged 
molecule potential energy, u. 





300 K, 


0.768 f 


5 cm 








u (kJ mol ^) 




{rihb) 


fo 


/i 


/2 


/3 


u > -48 




1.68 





0.28 


0.63 


0.05 


-48 < u > -103 




1.87 





0.19 


0.72 


0.08 


u < -103 




2.12 





0.04 


0.80 


0.16 


ensemble average 




1.86 


0.017 


0.19 


0.71 


0.082 






170 K, 


0.768 ^ 


; cm 








u (kJ mol ^) 




{nhb) 


fo 


/i 


/2 


/a 


u > -54 




1.91 





0.105 


0.88 





-54 < u > -108 




1.98 





0.05 


0.92 


0.028 


u < -108 




2.07 





0.0056 


0.92 


0.07 


ensemble average 




1.99 


0.00013 


0.04 


0.92 


0.037 



18 



Figure Captions 

1. Mean square displacements (MSDs) for different temperatures along (a) 0.768 g cm~^ and (b) 
0.878 g cm~'^ isochores. Note the logarithmic scale along both the axes of the plot. 

2. Static distribution, P{u), of tagged potential energy, u, of methanol, along (a) 0.768 g cm~^ and 
(b) 0.878 g cm~^ isochores at different temperatures. 

3. Contributions of reciprocal space, screened Coulomb and van der Waals interactions to the static 
distribution of tagged particle potential energy of methanol at 300 K and 0.768 g cm~^. P{urec) 
has been scaled by 1/12 while P{u^dw) has been scaled by a factor of half. 

4. Comparison of power spectra associated with temporal fluctuations in tagged particle potential 
energy, S„(f), between methanol (HI potential) and water (SPC/E potential) at 300K. The 
densities of methanol and water were taken as 0.768 g cm~^ and 1 g cm~^, corresponding to 
experimental densities at 1 atmosphere pressure. 

5. Power spectra of tagged particle potential energy fluctuations, 5„(/), of methanol, at different 
temperatures along (a) 0.768 g cm^^ and (b) 0.878 g cm^'^ isochores. Relative magnitudes of 
power spectra for different temperatures were scaled for clarity of presentation. 

6. Power spectra of tagged potential energy, Su{f ) of methanol fitted using the parameteric model 
defined in Section 2.3 along the 0.768 g cm~^ isochorc at (a) 300K and (b) 155 K. The raw power 
spectral data are shown with points and the overall fit is shown in solid lines. 

7. Effect of modifying the partial charge distribution of the HI model potential on: (a) Power 
spectra of fluctuations of tagged particle potential energies, Su{f)', (b) Static distributions of 
tagged particle potential energies at 300 K and 0.768 g cm~^. Relative magnitudes of power 

spectra, <S'„(/), are scaled for clarity. 

8. Arrhenius plot of self diffusion coefficient of oxygen along 0.768 and 0.878 g cm~^ isochores. 

9. Correlation plot between the self-diffusion coefficient of oxygen and the scaling exponent, of 
the multiple time-scale region of the Su{f) power spectra in methanol. Diffusivities are in units 
of cm^ sec~^. 

10. Correlation plot between the self-diffusion coefficient of oxygen in methanol and water, and silicon 
ion in silica with the scaling exponent, of the multiple time-scale region of the Suif) power 
spectra. Dimensionless diffusivities, D* = Dp^/^ / ^Jk^Tjm where p and m correspond to the 
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number density and masses of the formula units have been used to facihtate comparison between 
the different systems. 
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Figure 1: 
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Figure 2: 
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Figure 3: 
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Figure 8: 
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Figure 9: 
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Figure 10: 
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